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Abstract 

We study the renormalization group equations of Ma's scotogenic model, which 
generates an active neutrino mass at 1-loop level. In addition to other benefits, the 
main advantage of the mechanism exploited in this model is to lead to a natural loop- 
suppression of the neutrino mass, and therefore to an explanation for its smallness. 
However, since the structure of the neutrino mass matrix is altered compared to 
the ordinary type I seesaw case, the corresponding running is altered as well. We 
have derived the full set of renormalization group equations for the scotogenic model 
which, to our knowledge, had not been presented previously in the literature. This 
set of equations reflects some interesting structural properties of the model, and it 
is an illustrative example for how the running of neutrino parameters in radiative 
models is modified compared to models with tree-level mass generation. We also 
study a simplified numerical example to illustrate some general tendencies of the 
running. Interestingly, the structure of the RGEs can be exploited such that a 
bimaximal leptonic mixing pattern at the high-energy scale is translated into a 
valid mixing pattern at low energies, featuring a large value of 6*13. This suggests 
very interesting connections to flavour symmetries. 
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1 Introduction 



The leptonic sector of the Standard Model (SM) of particle physics has undergone several 
modifications within the last decades. Not only have we learned from oscillation experi- 
ments that neutrinos must have a non-zero mass, but we have meanwhile also measured 
the corresponding mixing angles [TJ [2] . This includes in particular the very recent spec- 
tacular measurements of the previously unknown mixing angle 6*13 by the Daya Bay [3] 
and RENO $\ collaborations, after hints from MINOS [5], T2K [6], and Double Chooz 
What is still missing, however, is the knowledge of the absolute neutrino mass scale: Al- 
though we have limits from kinematical measurements such as Troitsk [8] and Mainz [9], 
from experiments on neutrino- less double (3 decay [TOT [TT] . and from cosmology [12], the 
only decent statement we can make is that the true neutrino mass scale is somewhere 
below 1 eV. 

Nevertheless, generating a mass scale < 1 eV in a model is far from trivial: In the 
classic version of the SM, neutrinos are strictly massless, but even if we extend the SM 
by right-handed neutrinos to generate Dirac neutrino mass terms, the tininess of the 
neutrino mass scale compared to the scale of electroweak symmetry breaking (EWSB) 
looks highly unnatural. Probably the most frequently discussed solution to this problem 
is to introduce a Majorana mass term for the right-handed neutrino states and to exploit 
the potential large size of the corresponding masses through the famous seesaw (type I) 
mechanism jT3j [TTJ [151 HI Eh as depicted on the left panel of Fig. [TJ Seesaw-type 
settings have also been investigated in the context of renormalization group equations 
(RGEs) [TBI CS1 I2D1 ED [221 122], in order to get a more detailed picture of the change 
of neutrino masses and mixing parameters with the energy scale, usually referred to as 
RG-running or simply running. 

An alternative and very attractive way to explain the smallness of neutrino masses 
is to assume that they are strictly zero at tree-level, but receive non-zero higher order 
corrections. Maybe the most simple such framework can be found in Ma's scotogenic 
model [21] (for simplicity called "Ma- model" in this paper), in which the introduction 
of only a few new fields and one additional symmetry compared to the SM leads to a 
non-zero neutrino mass at 1-loop level, cf. right panel of Fig. [TJ In fact, this diagram 
exhibits a structure that is somewhat reminiscent of a seesaw-type situation, which is 
why sometimes the corresponding suppression is dubbed as radiative seesaw. 

The main idea that led to this paper is to investigate the running of neutrino param- 
eters in models with a radiative neutrino mass, as the rich structure of the corresponding 
diagrams is likely to introduce interesting and yet unknown features in the running, which 
are not present in tree-level realizations of the seesaw mechanism. To our knowledge, there 
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Figure 1: Tree-level and radiative seesaw mechanisms. 



exists no such study in the literature with the focus put on the neutrino sector in radiative 
models, and we aim to start this enterprise by a study devoted to the RGEs of the Ma- 
model. Naturally, this could be extended to other radiative models for neutrino masses, 
such as the Zee-Babu model [251 [26] or the Aoki-Kanemura-Seto model [271 EE]- In par- 
ticular the interplay between the scalar and the lepton sectors has the potential to reveal 
interesting new effects, as we will already see in this study. 

However, we want to stress that several studies are already available which investigate 
e.g. limiting cases of our framework or subsets (or generalizations of subsets) of certain 
sectors of the Ma-model. A particular example for such a case would be the investigations 
of the RGEs of a general Two Higgs Doublet Model (THDM). Whenever applicable in 
this paper, we will refer to the corresponding works treating these related frameworks. 

This paper is organized as follows: In Sec. |2j we review Ma's scotogenic model and 
discuss the different effective theories arising when subsequently integrating out the heavy 
neutrino fields. Next, in Sec. El we discuss in detail the matching conditions at the 
boundaries between the respective theories, which in our case have to be consistently 
imposed at 1-loop level. Our main results, the explicit RGEs at 1-loop level are presented 
in Sec. HI After that, we present a numerical exemplifying study (in a slightly simplified 
framework) in Sec. El in order to illustrate how to use our results. We finally conclude in 
Sec. El 



2 Ma's scotogenic model 

The so-called scotogenic model has been discussed by Ma [23] , and in the following we will 
therefore call it Ma-model for simplicity. In this section, we will first review this model, 
and then discuss some of its low-energy limits, which we will also use in our calculations 
later on. 
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2.1 The Ma- model as full theory 

The Ma-model is a very economic extension of the SM, which however exhibits a vari- 
ety of very interesting features. Some of its nicest aspects are its possibilities for Dark 
Matter (DM) candidates, the radiative generation of neutrino masses, and the collider 
phenomenology of relatively light right-handed neutrinos. Compared to the SM, the 
model contains three new crucial ingredients: 

1. In addition to the SM Higgs there is a second scalar doublet rj ~ (2, 1/2), which 
has the same quantum numbers under SU(2)l x U(1)y as the ordinary Higgs. 

2. Furthermore, there are three right-handed neutrinos Nk ~ (1 ; 0), which have a 
Majorana mass term, £ = — ^M^N^N^ 

3. Finally, there is an exact (discrete) Z 2 parity, under which all SM-fields transform 
trivially (i.e., they have ^-charge +1), while all new fields, rj and Nk, are odd (i.e., 
they have ^-charge —1). Because of the non-trivial charge under this symmetry, 
the new doublet rj is often dubbed as inert or dark scalar doublet. 

These few assumptions lead to a variety of very interesting consequences, some of 
which we will list in the following: 

• Due to the Z 2 -symmetry being unbroken, the lightest particle charged under this 
symmetry will be absolutely stable. If this particle is electrically neutral, i.e., it is 
one of the Nk or one of the neutral components of the scalar field rj, it would be a 
good candidate particle for DM [291 EQj EJ Ml 1321 E3J- 

• Since the standard neutrino Yukawa coupling, £Y U k,v = — LY U §N R + h.c. where 
$ = i<T2$*, is forbidden by the exact Z 2 symmetry, this term cannot lead to a tree- 
level neutrino mass. There is an alternative Yukawa coupling involving the new 
scalar, £ Yuki , = ~Lh u fjNR + h.c, but also this term does not generate a neutrino 
mass at tree-level. However, it does lead to a non- vanishing neutrino mass at 1-loop 
level, cf. right panel of Fig. [TJ The resulting formula for the neutrino mass matrix 
is: 



16n 2 



hikhjk njT m\ ^ m| _ mj ^ mj 
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This expression involves the scalar masses m# (mj) of the real (imaginary) parts 
of the electrically neutral component 77° of the new scalar. These masses originate 



1 Note that we have already used the freedom of rotating the basis such that this Majorana mass term 
is diagonal. 



3 



from the scalar potential: 



1 2 1 

^scalar = m\& <$> + m^T] + -Ai ($ t $) + -A 2 (^) ' 



+ A 3 ($ f $) (77*77) + A 4 ($ t r/) (77+$) + -A 5 ($+77) 2 + h.c. 



(2) 



where all parameters can be taken to be real, but we need to choose m\ < and 



m\ > in order to have a VEV of ($°) 



174 GeV only, while keeping (77°) = 0. 



This choice leads to the following scalar masses (just as for a general THDM 
with A 6j7 = 0), where h is the SM-like Higgs particle: 

m\ = 2\ lV \ 

m R = ^(Vr) = m l + (-^3 + A 4 + A 5 )t; 2 , 
m] = m 2 (r]°j) = m\ + (A3 + A 4 - X 5 )v 2 , 
m 2 (r] ± ) = ml + X3V 2 . 



(3) 



The formula for the light neutrino masses, Eq. (JT]), exhibits a behaviour similar to 
that of a mass matrix resulting from a seesaw mechanism, since the Dirac Yukawa 
couplings are suppressed by some high mass scale. Hence this mechanism is often 
referred to as radiative seesaw. However, due to the additional loop suppression 
factor (167T 2 )" 1 , and due to a hidden proportionality to the small coupling A5, only 
part of the suppression has to come from the masses of the heavy particles involved. 
Accordingly, the masses of the right-handed neutrinos can have relatively low values, 
i.e. around the TeV-scale, similar to the masses of the new scalar, thereby leading 
to potentially interesting collider signatures [521 EHJ E3 155] . 

• Due to the extended scalar sector, this model typically leads to potentially observ- 
ables rates of lepton flavour violating (LFV) processes [2SJ E2] • In turn, when trying 
to explain a more detailed structure in the flavour sector, exactly these LFV pro- 
cesses tend to spoil the application of flavour symmetries to predict observables like 
mixing angles |39j . 

• On the other hand the extension of the Ma-model by an additional Left-Right 
symmetry can yield a completely new mechanism to explain large but not necessarily 
maximal mixing angles, the so-called radiative transmission of hierarchies [40]. This 
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mechanism could be an alternative to flavour symmetries and in that way resolve 
the problem mentioned in the previous point. 

As we can see, there are indeed many interesting aspects of this model discussed in 
the literature. What is, however, to our knowledge lacking up to now is a thorough 
treatment of the running of the neutrino mass matrix in the Ma-model. Since running 
is unavoidable, we should carefully check in how far the behaviour changes in radiative 
models, like the one under consideration, compared to models with tree-level neutrino 
masses. We could expect at least some qualitative changes, due to different diagrams 
and structures appearing in the calculation. Furthermore, since the neutrino mass is zero 
at tree-level, its 1-loop correction is actually the decisive term, and the running of this 
term could be potentially large. It is therefore important to carefully check if consistency 
of this model with low-energy data arises naturally, or if it maybe could be under some 
pressure due to the running effects being too strong. 

2.2 The effective field theories appearing 

In particle physics, it is well-known that a separation of scales translates into a decou- 
pling of the heavy degrees of freedom at low energies, which goes under the name of 
the Appelquist-Carazzone theorem [U], leading directly to so-called effective field theories 
(EFTs) [12]. In practice, this mans that we can encode the effects of new physics at high 
energies into effective operators, which contain only fields of the low-energy theory, but 
which have a mass dimension d larger than four. In an EFT, we solely take into account 
the terms in the Lagrangian related to fields which are relevant at the corresponding 
energy scale, i.e., we drop the terms involving the momenta of very heavy fields. It is this 
procedure that leads to the higher- dimensional terms. Note that there is an upper bound 
on the validity of any effective field theory, usually denoted as the cutoff A, inverse powers 
of which suppress the effective operators. Actually, EFTs are just a way to "parametrize 
our lack of knowledge" of the underlying low-distance physics. 

We now take on a top-down approach, in which we extract the EFT from a more fun- 
damental high-energy theory. This means that we start with a completely renormalizable 
high-energy theory, run it down to lower energies by gradually decoupling the heavy fields, 
and eventually obtain the low-energy footprints of the high-energy regime of this theory. 
In order to do this we make use of the corresponding renormalization group equations. 
In the case of the Ma-Model, we will consider the decoupling of the heavy neutrinos - 
they will be integrated out - and we will investigate which terms appear in the result- 
ing low-energy effective Lagrangians. Note that we have assumed the inert scalar to be 
lighter than the right-handed neutrinos for the sake of an example. Then, in a manner 
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very similar to a type I seesaw mechanism, integrating out the heavy neutrinos will lead 
to d = 5 operators coupling two lepton doublets to two scalar doublets. 

To have a consistent way to describe the running of the neutrino masses, we hence 
make use of the so-called Weinberg operator [43J , 

£f 5 = \* 8 f Wc^J f Lb e ba <p a +h.c. , (4) 
4 v ' 

=O k : effective operator 

where e^- is the two-dimensional antisymmetric tensor. Indeed, the 1-loop diagram in 
Fig. [1] contributes to the Weinberg operator, and it corresponds to the type T-3 contri- 
bution in the general classification in Ref. jH]. With only the SM particle content, the 
Weinberg operator is the only d = 5 operator possible, and it will give masses to neutrinos 
once the Higgs field obtains a VEV. However, in our case we have an additional scalar 
doublet r], which allows for more Weinberg-like operators: 

r,s=l 

where (jP^ = $ and (f)( 2 > = rj. These four potential operators are depicted in Fig. [21 
However, it is clear that the Z 2 symmetry of the Ma-model forbids two of these oper- 
ators, and Ok , and thus the final remaining operators are 0i n) and C>! 22) . For 
convenience these operators will be written as and , but one has to keep in mind 
that this is not purely correct, as the k's are actually just the coefficients of the d = 5 
operators, and not the operators themselves. Note that will not directly contribute 
to the neutrino mass, as the inert scalar does not obtain a VEV. However, it will mix with 
the SM Weinberg operator at 1-loop level, cf. Fig. [3], making it necessary when attempting 
to impose a consistent matching at 1-loop level (cf. Ref. [45J). As a result, the RGEs of 
these two operators are coupled. 

Now we know that our low-energy effective theory contains two effective operators 
(which are mixed at 1-loop level) contributing to the mass of the neutrinos. We can then 
renormalize the effective theory and evolve our set of parameters (Yukawa couplings, scalar 
couplings, scalar masses, effective operators) to the scale at which this theory breaks down, 
typically the mass of one of the heavy neutrinos. However, there is no reason why the three 
neutrino singlets would be degenerate in mass. Then, instead of integrating out all the 
heavy neutrinos at once to go directly from the full to the lowest-energy effective theory, it 
would seem more natural to include a certain hierarchy among the heavy neutrino masses 
and to integrate out the corresponding fields one by one. Consequently, we will obtain a 



= E +h-c , (5) 

=0k .effective operators 
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Figure 2: The potential Weinberg-like operators in a general THDM. In the Ma-model, 
the two operators and are forbidden by the unbroken Z 2 symmetry. 
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Figure 3: The 1-loop mixing of the two Weinberg-like operators in the Ma-model. 



set of three different effective theories characterized by their own respective scales, i.e., 
the three masses of the heavy neutrinos. These Majorana masses are unrelated to the 
EWSB scale and they will be considered to be of the order of several TeV. Later, the 
three effective theories will be denoted by EF1, EF2, and EF3, while the full theory will 
be denoted by FT. 



3 Matching the theories 

Let us summarize what is needed to be done to get the global running of the parameters 
over the set of theories: 

1. We start at a scale above the masses of all the particles, where the renormalizable 
theory (i.e., the FT) is valid. 

2. The next step is to evolve this theory down to lower scales. As long as the energy 
scale under consideration lies above all heavy neutrino masses, this evolution is 
described by the RGEs of the FT. 
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3. When the energy scale considered goes below the mass of one of the neutrino 
singlets in the theory, we must change the theory to the new theory ETA: without 
the corresponding field N^, thereby introducing non-renormalizable interactions via 
effective operators. 

4. Both the changes in the existing parameters and in the coefficients of the new 
interactions are computed by matching the physics directly above and below the 
boundary between the two theories. This must be done at each mass scale boundary. 

To achieve the first three points, it is just necessary to renormalize each theory, tak- 
ing into account the changes in the particle content and in the effective operators. We 
will end up with parameters whose running is described by piecewise continuous func- 
tions. To avoid a mismatch at the junctions between theories, we need the fourth point: 
Indeed, as often in physics, the use of highly simplified models brings some unphysical 
behaviour. Here it is obvious that the heavy neutrinos do not decouple abruptly, but in- 
stead their influence weakens slowly until it becomes insignificant. The exact behaviour is 
complex, but we can fix the potential mismatch of the running between different theories 
by imposing so-called matching conditions. This procedure has the advantage of being 
really straightforward, and it consists only in evaluating both theories at the matching 
scale, i.e., the energy scale at the boundary. This procedure is standard in the study 
of the running of the parameters when considering different effective theories (see, e.g., 
[T8~| H~9| |23| |4"5| H6]). However, it is usually performed at tree- level, whereas we must 
consider it at 1-loop level to obtain a consistent answer, due to the 1-loop mixing be- 
tween fi^ 11 ) and kS 22 \ The corresponding diagrams for the different effective theories are 
depicted in Tab. [TJ Power counting indicates that they do not all have the same degrees 
of divergence. Indeed, diagrams in ETs have typically higher degrees of UV divergence, 
as they contain fewer propagators. For example, in the case of k^ u \ the second diagram 
in ET3 is logarithmically divergent, while the corresponding diagram in the FT is finite. 
This is not an obstacle, but we will have to regulate each diagram using dimensional reg- 
ularization. Once that is done, the infinite part will be skipped, as it can be compensated 
by an adequate counterterm in the Lagrangian. Then, only the finite scale-dependent part 
will matter for the matching. Note that, in our case, we have one more interesting sim- 
plification: Due to the 1-loop diagrams being proportional to some parameters Aj from 
the Higgs potential, which is not the case for the tree-level diagrams, we can separate 
the tree- and 1-loop level contributions and match them separately rather than matching 
their sums. The matching scheme for the tree-level diagrams is depicted in Tab. [2] for 
illustration. 
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Table 1: Overview of the diagrams to consider for the matching in the different theories 
(FT & ET3). Effective operators are drawn in a space-saving version. 



We have performed the calculations of all these diagrams, and of all other diagrams 
relevant for the running. This includes several running quantities such as gauge couplings 
or parameters from the Higgs potential, above the energy scale of EWSB, where only the 
scalar doublets and the neutrino singlets have mass. Note that, at the lowest energies, 
this actually introduces a slight inconsistency: EWSB would set in, all particles would 
become massive, and in particular we might have to integrate out some of the inert scalars. 
However, that would happen at a scale very close to Mz, and unless the threshold effects 
are extremely strong the resulting effect would be negligible. We will leave the detailed 
investigation of this parameter region for future studies, as the main results presented in 
this paper are the RGEs and the discussion of one very illustrative case. 
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Table 1: (continued) Overview of the diagrams to consider for the matching in the different 
theories (ET2 & ET1). Effective operators are drawn in a space-saving version. 



We now impose the matching conditions. This is a systematic procedure which is 
done order by order in the loop-expansion. When two theories are compared at a given 
loop order, the lower order results have to be included in the matching. In our special 
case, however, they happen to decouple. Furthermore, despite of having a stack of three 
effective theories, going from one to another is each time a similar process, namely the 
decoupling of exactly one heavy neutrino (unless two or more of them were degenerate 
in mass). Therefore, the matching conditions will be analogous and we can write down 
a general expression, valid whichever the ET relevant for the energy scale is considered. 
For simplicity, we adopt the convenient notations from Ref. [IS] : 

• /X r ) — y ( K,(rr) is the effective operator in the n-th effective theory, i.e., the op- 
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Table 2: The matching of the tree-level diagrams. None of these operators can directly 
lead to a light neutrino mass, but they will contribute at 1-loop level. The abbreviation 
"MC" stands for matching condition. 
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erator making up for the low-energy influence of the heavy neutrinos TV*, (with 
k = n, . . . , 3). 

Mtv — >■ Mtv is the mass matrix of the remaining heavy neutrino singlets in ETn. 
In the FT this matrix is of dimension 3x3, but in ETn the states Nf, with k > n 
have been decoupled, and their masses are not defined anymore. As a result the 
heavy neutrino Majorana mass matrix in ETn is an (n — 1) x (n — 1) matrix. We 

(n + l) 

denote by M n the largest eigenvalue of M . 

(n) 

h v — > h v is the neutrino Yukawa matrix taking in account that we have integrated 
out (4 — n) generations of neutrino singlets. It implies that it is a non-square matrix 
coupling of (n — 1) neutrino singlets to 3 lepton doublets. Taking rows corresponding 
to neutrino singlets and columns corresponding to lepton doublets, the coupling 
technically looks like an (n — 1) x 3 submatrix of the "complete" Yukawa matrix: 



h u 



(Ml,2 (Ml,3 \ 

(hv)n-l,\ {hv)n-l,2 {h u )n-l,2, 
••• 



(n) 

\ y h v 



Note that with these notations, we have: 

(1) (rr) (" >3 > n ( u u ( u n fa\ 

K( rr ) = K J , K(f. r ) = U, h v — h v , h v = U. (oj 

In practice we merely match the values of the diagrams, which change when going from one 
theory to another. Considering the scheme from above and performing the calculation, 
one obtains the following matching conditions (with n = 1, 2, 3): 

1. For k^: 

MCn(/j, = M n ) — > K (11) , y = K (11) , y + —— Ah v ) in M n ( /i „) fti . (7) 
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2. For k^: 



MCn (ji = M n ) 



(n) 



(n+1) 



K (22),ij — K (22),ij T 



2 + 



f{M n ,m 2 ) 
g(M n ,m 2 ) 



(1+1) 



(n+1) 



ln(^ 



7S and 



Here, we have used the loop functions and f(x, y) = ^_ y ^ - 

In ^^2) ) which arise from the computations of the corresponding diagrams [IT]. In general, 
H denotes the renormalization scale. 

There is one more subtlety involved: The matching conditions, as presented here, are 
valid only in the basis where the heavy neutrino Majorana mass matrix is diagonal, which 
is the setup we imposed at the input scale. However, the running does not care for our 
choice of basis and it will introduce some off-diagonal coefficients in our initially diagonal 
matrices [19]. Therefore, each time when we integrate out one heavy neutrino at the 
corresponding threshold, the Majorana mass matrix will have to be diagonalized first, so 
that we can use the derived expressions. This procedure imposes a redefinition of the 
Yukawa matrix that can be written as follows: 



(n) 

M 



(n) 



VM^V 1 



(») 
h v 



V T h v . 



(n) 



(9) 



Here, V must be a unitary matrix since M is Hermitian. We have to perform this "re- 
diagonalization" at each matching scale. 

Finally, we can write down the most general expression for the light neutrino mass 
matrix, arising from all contributions (only at 1-loop in this case) to the operator k^ 11 ': 



A, 



16tt 5 



g(li,m 2 



(«) («) (n) 

« (22) + h u M^ 1 f{M n ,m 2 )h l 



(10) 



Note that the structure of this matrix is actually very similar to the one in an ordi- 
nary seesaw type I framework (or in the corresponding effective theories, respectively), 
but corrected by some loop functions. To determine the leptonic mixing, it is enough 
to diagonalize this light neutrino mass matrix with a unitary matrix t/^/i), as well as 
the charged lepton mass matrix with a unitary matrix U e (fi). Note that these matri- 
ces, naturally, depend on the energy scale /i under consideration, which will in the end 
translate into the running of the neutrino mixing parameters. Hence, the Pontecorvo- 
Maki-Nagakawa-Sakata (PMNS) leptonic mixing matrix is directly given by: 



PMNS 



11] 
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We adopt for U = ?7pmns the following parametrization: 



C12C13 



S12C13 



U 



~Sl2 c 23 — Ci2S23Sl3e* 5 C12C23 — Si2S23 s 13 el<5 S23 c 13 



S12C23 - Ci 2 S23Si 3 e 



-C12C23 - s 12 s 2 3S 13 e t5 C 23 Sl3 




(12) 



where Sij = sin and c,j = cos 6^-. The mixing parameters can then be extracted from 
the PMNS-matrix by 

6»i3 = arcsin(|E/i 3 |) , 







arctan 



12 



\U 



12 



t/ll 



arctan 



7 23 



\u. 



23 | 



33 | 



if U n + 
else 

if U 33 ± 
else 



-arg 



2 

Ci2C 13 C23Si3 



+ C12C23S13 / (S12S 



'23y 



(13a) 
(13b) 

(13c) 
(13d) 



Now that we have discussed the principle ideas, we can move on to our results, by first pre- 
senting the analytical forms of the RGEs and then investigating their numerical solutions 
in a simplified example scenario. 



4 One-loop Renormalization Group Equations 

The main result of this work is the set of 1-loop RGEs for Ma's scotogenic neutrino mass 
model, to be presented in this section. The presentation of this full set of RGEs (with 
the only exception of the effective theory with the scalars being integrated out) is, to 
our knowledge, a completely new result. However, of course certain limiting cases have 
already been discussed in the literature. 

Renormalization group equations have been extensively studied, e.g., in the contexts 
of the SM and the MSSM. The references to these are numerous and sometimes messy, so 
we will refer to the reader only to the very nice Refs. [211 SSI ED] for the presentation of 
an exhaustive set of RGEs. However, in the theoretical context of the model investigated 
here, these equations are not adequate, and the best way to find consistent results is to 
consider the derivation of the RGEs in a general quantum field theory, as performed in 
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Refs. [5H Ml E3], or in a general THDM [51]. 

More specifically, type I seesaw models have been discussed in Refs. [HI [191 [20] . Tree- 
level extensions of the seesaw mechanism have been discussed in Refs. [22j [23j H6j HH] , and 
the case of Two- (or Multi-) Higgs Doublet Models has been presented in Refs. [TBI 155] . 
Stability of texture zeros has been investigated in Refs. [2U [56]. Studies of the Weinberg 
operator (or Weinberg-like operators) have been performed, e.g., in the context of the 
SM [57] or in a general THDM [15]. Finally, extremely useful software packages have 
been developed in connection with Ref. [H]: REAP ( "Renormalization Group Evolution 
of Angles and Phases") [58] and MPT ("Mixing Parameter Tools") [59], to be found 
online. Whenever applicable, we have checked the consistency of our results with these 
references, and we have always been able to reach agreement. 

Hereafter we present the set of RGEs corresponding of the Ma-model, where through- 
out the paper we will use the abbreviation 

V = IQtt 2 ^, (14) 
d/i 

with /i being the renormalization scale. Then, the following parameters will be involved 
in the running, and their running (in the FT) is determined most easily in the order 
indicated here: 

1. The gauge couplings gi, where i = 1, 2, 3. 

2. The SM-like Yukawa coupling matrices Y x , with x = u,d, e, for the up-quarks, down- 
quarks, and charged leptons, as well as the neutrino Yukawa coupling matrix in the 
Ma- model, h v . 

3. The heavy neutrino Majorana mass matrix and the five real scalar couplings 
Aj, with i = 1,2, 3, 4, 5, which appear in the Higgs potential. 

4. The mass parameters mi and m 2 of the SM-like Higgs and of the inert scalar, 
respectively. 

In addition to these quantities, we also have the effective operators «(n) and ( AC(22) 
running in the different effective theories ETn, whose RGEs are coupled, cf. Sec. 12.21 All 
these RGEs have been calculated using dimensional regularization, while adopting the 
minimal subtraction (MS) scheme. We will now present the 1-loop RGEs and discuss 
the different sets of equations one by one. While in this paper we present the resulting 
equations, more details on their derivation and on the actual calculation can be found in 
the corresponding thesis, Ref. [47] . 
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Figure 4: The runnings of the gauge couplings and of the top Yukawa coupling in the 
SM and in the Ma-model. 



4.1 The gauge couplings 

The RGEs for the gauge couplings are given by: 

Vg 2 ~- 
Vg 3 = -7gl 



3gl 



(15a) 
(15b) 
(15c) 



Compared to the SM, these equations are only slightly altered because they do not depend 
on the Yukawa and scalar couplings but only on the number of fermion generations and 
the number of scalar doublets in the theory [60]. Note that the differential equations fl!5a[) 
to f)15cl) can all be easily solved analytically, leading to 



9i(fk>) 



lMVo)ln 



(16) 



where (61,62,^3) = (7,-3,-7) and fiQ denotes the reference scale. The running of the 
gauge couplings is depicted on the left panel of Fig. HJ In particular the running of g 2 
shows a considerable difference in the Ma-model compared to the SM, due to the presence 
of an additional scalar doublet. Note that we have chosen the input values for the gauge 
couplings such that the measured values are correctly reproduced at the Z-pole, i.e., at 
the scale \i = Mz- 
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4.2 The Yukawa couplings 

The RGEs for the Yukawa coupling matrices are given by: 

VY U = Y u { \YlY u - h}Y d + T- l lg\ - \g\ - %g\ } , (17a) 
VY d = Y d S^YjY, - 3 -Y}Y u + T- - \g\ - 8, 3 2 } , (17b) 
VY e = Y e !^YX + \hlh u + T — ^gl - 9 -gj} , (17c) 

Vh u = h v \^-hlh v + l -Y}Y e + T V - \gl - ^ 2 2 | , (17d) 

where T = Tr (y}Y c + 3YjY u + 3Y}Y d ^j and T„ = Tr (hlK). Note that, compared to 
the case of the SM supplemented by a Dirac neutrino Yukawa coupling matrix Y u , the 
neutrino Yukawa coupling matrix h u in the Ma- model only appears in Eqs. f)17cp and (I17dj) . 
However, numerically Eqs. (I17ap and ( 11Tb[) will not really differ from their equivalents in 



the SM with Dirac neutrino masses, due to the extremely small entries in the Yukawa 
matrix Y v . The running of the top Yukawa coupling, in the approximation of this being 
the only non-zero entry of Y u , is depicted on the right panel of Fig. HI As to be expected, 
the top-running is dominated by the diagrams involving strongly interacting particles, 
and hence it is not really different in the Ma-model compared to the SM. Note that the 
set of equations (I17ap to (I17dj) can be solved independently of all other parameters but 
the gauge couplings. Hence, by inserting the solutions of the equations from Sec. 14. 1[ the 
complete evolution of the Yukawa couplings can be studied. 

4.3 The right-handed neutrino mass matrix 

The RGE required for the running of the mass matrix of the neutrino singlets can be 
found, e.g., in Refs. [T9| l6Tj : 



VM N = {M N (Khlf + (h u hl) M^v} . (18) 

This equation is quasi identical to its seesaw type I equivalent, cf. Ref. [20]. Note that the 
only "external" input in this RGE is the Yukawa coupling matrix of the neutrinos, h u . 
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Hence, once the solution to Eq. (117d|) is known, we can immediately obtain the numerical 
evolution of M N . 



4.4 The scalar sector 

The RGEs of the scalar couplings have been calculated for the SM and its extended 
version including heavy neutrinos [TSJ [2TJ ED]- I n the context of a theory with more 
than two scalar doublets there are less studied, but some references give very general 
results [54] . However, in our specific case involving an additional scalar doublet as well as 
three neutrino singlets and an exact Z 2 symmetry, these results are not adequate. Thus, 
we had to adapt the RGEs to our case, leading to: 



VX 1 = 12X1 + 4A^ + 4A 3 A 4 + 2X 2 4 + 2X 2 5 + - (g* + 3g\ + 2g\g 2 2 ) 



(19a) 



-3Ai 



g[ + 3g\ - -T 



4T 4 



VX 2 = 12X1 + 4A§ + 4A 3 A 4 + 2X 2 4 + 2Xj + - (g{ + 3g\ + 2g\g 2 2 ) 



-3A S 



g\ + 3^ - ~T V 



-4T, 



An ■ 



(19b) 



VX, 



VX A 



(Ax + A 2 ) (6A 3 + 2A 4 ) + 4A* + 2X\ + 2X\ + - {g\ + 3g\ - 2g\g 2 2 ) 



-3A, 



g\ + ?>gl " 77 ( T - + T ) 



AT 



2 (Ai + A 2 ) A 4 + 8A 3 A 4 + 4A* + %X\ + 3g 2 g 2 2 

2 



-3A 4 



g{ + m - o % + T ) 



4T„ 



VX 5 = 2 (Ai + A 2 ) A 5 + 8A 3 A 5 + 12A 4 A 5 - 3A E 



gl + 3g 2 2 --{T v + T) 



(19c) 



(19d) 



(19e) 



where T 4 = Tr {y^X + 3YjY u YjY u + 3^^]^}, T 4v = Tr {hlKhtK}, and T„ e = 
Tr {h^JiyY^Y^ . These are the most demanding equations in practice, since one needs the 
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knowledge about the evolutions of all gauge and Yukawa couplings as input. Nevertheless, 
after having completed the program from Sees. I4.11 and I4.2[ it is no problem to also obtain 
the evolutions of the 4-scalar couplings Aj. 

Note that the corrections to all scalar couplings except for A5 are additive, while A5 
only receives multiplicative corrections. This does not happen by accident, but rather 
there is a very deep reason for that: The Ma-model intrinsically involves the breaking 
of lepton number either by the heavy neutrino Majorana mass term or directly by the 
As-term, cf. Sec. 12. 1[ depending on which field, Nk or 77, lepton number is assigned to. 
However, even in the case of lepton number violation coming in directly by the heavy 
neutrino mass term, this breaking is not easily transmitted to the active neutrino sector. 
The decisive point is that the Dirac neutrino Yukawa coupling, C Yuku = —Lh v fjNn + h.c, 
is lepton number conserving. However, the diagram giving mass to the light neutrinos, 
cf. right panel of Fig. [TJ and the corresponding neutrino mass matrix, Eq. (fIU|) . are 
proportional to A5, and they will transmit the violation of lepton number from the heavy 
to the light neutrino sector. Hence, in the case of A5 = 0, either the lepton number is 
exactly conserved (if n carries lepton number) or violation is present in the model but will 
not translate into the light neutrino sector (if the carry lepton number). In any case, 
the form of Eq. ( I19ef) is a protective limit of lepton number symmetry: If A5 = holds at 
the input scale, it will always remain true, at any energy. This enforces all corrections to 
A5 = to vanish identically, which is only possible if they are proportional to some power 
of A5. In addition, supposed that we start with very small couplings (close to zero) at 
the high-energy scale, one could hope for the parameters Xi,2,3,4 to naturally grow to the 
relatively large values needed, while A5 remains very small. However, as we will see, at 
least for A2 such a magical behaviour is destroyed by the requirement of keeping vacuum 
stability. On the other hand, it is still true for Ai^s. 

We have also calculated the RGEs for the scalar masses in the Higgs potential: 

Vm\ = 6Aim? + (4A 3 + 2A 4 ) m\ + m\ 

Vm\ = 6\ 2 ml + (4A 3 + 2A 4 ) m\ + m\ 

Note that the neutrino Yukawa coupling matrix only contributes to m\, while all other 
Yukawa couplings contribute only to m\. Furthermore, we also have to input the evo- 
lutions of Ai j2 ,3,4 and of the gauge couplings <7i ;2 into the equations for the scalar mass 
parameters. 



2T 



3$ 



(20a) 



2T„ - - (g\ + 3g 2 2 ) 



(20b) 
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4.5 Running in the effective theories 

Finally, we have computed the RGEs in the different effective theories. As explained, 
when integrating out the heavy parameters of the theory (chosen to be the heavy neutrino 
singlets), we switch from the full theory to the tower of effective theories, and the set of 
RGEs is altered compared to the FT. Hence, we need RGEs for the d = 5 effective 
operators which generate a neutrino mass at low energies. 

In the context of the THDM, such an operator is potentially split into four parts [45l 
However, as explained in Sec. 12.2} the Z 2 symmetry allows only two terms which we 
renormalize. The /3-functions of the remaining two operators arej^l 

Vk^ = -~ {^ U \YX) + (l?y e ) T « (11) } + 2T K ( n ) + 2A 1 k( 11 > + 2\ 5 k^ - 3<7 2 V n ), 

(21a) 

Vk^ = I {^ 22 \Y}Y e ) + (y e ty e )V 22 )} + 2\ 2 k^ + 2A 5 « (11) - 3 9 2 2 k^. (21b) 

This allows to write the general form of the RGEs throughout the different effective 
theories. 

5 The numerical solution of the RGEs 

Finally, we present a numerical analysis of the equations presented in Sec. SJ Since the 
whole set of RGEs is rather complicated, and a detailed numerical study of many example 
setups is beyond the scope of this paper, we chose to investigate one approximate case to 
illustrate the most generic features of the RGEs of the Ma-model. Note that this case is 
not more than an example, and even if we have a fair agreement with experiments, we 
would still need a more complete investigation in order to make solid predictions. We 
propose such a detailed numerical investigation for future work. 
In our calculation, we make use of the following approximations: 

• The top-quark Yukawa coupling being much larger than the other Yukawa cou- 
plings, we will assume that down-quark and charged lepton Yukawa matrices can be 
neglected compared to the up-quark Yukawa matrix, of which only the top-quark 
component will be considered: Y& , Y e -C Y u ~ diag(0, 0, y t ). 

2 This result is in agreement with the more general case computed in Ref. [55| . However, if one wants 
to compare our results with the ones presented there, it is important to take into account the differences 
in the definitions of the coefficients in the scalar potential. 
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• The down-quark and charged lepton Yukawa matrices will be neglected compared 
to the gauge couplings: , Y e <C g iy where i = 1,2, 3. 

• We consider the charged lepton mass matrix to be diagonal. This simplification is 
motivated by the weak running of the charged lepton masses, and for the sake of an 
easy example one can accept the small errors introduced by it. 

In addition to these approximations, we also make certain simplifying assumptions: 

• As has been said previously, we assume that the inert scalar mass parameter (7722) 
is below the masses of the heavy neutrinos (Mi^), so that we can integrate these 
out while keeping the inert scalar as dynamical particle in the model. We take the 
physical mass to be m 200 GeV at the Z-pole, which is in agreement with the 
upper bound derived from LHC measurements [63]: m, q ± < 200 GeV. 

• We take the neutrino Yukawa coupling matrix h u to have real entries only, for 
simplicity. 

• We assume normal ordering for the light neutrino masses. 

• We assume that there is bimaximal mixing at the GUT (Grand Unified Theory) 
scale, M GUT = 10 16 GeV. 

Equipped with these approximations and assumptions, we will perform the numerical 
study using the top-down approach, i.e., we will run the parameters from high-energy to 
low-energy, using input values at the GUT scale. As outlined in Sec. HI the set of RGEs 
will be solved step by step starting with the simplest non-coupled ones. The evolution 
of the gauge couplings and of the top- Yukawa coupling had already been presented in 
Fig. HI These solutions are independent of the particular setting considered for the neu- 
trino sector, and can hence taken to be universally truejf] However, what we are mainly 
interested in is the neutrino sector and its interplay with the scalar sector, which will be 
investigated in the following. 

5.1 The evolution of the couplings 

Following our previous approximations, in the FT, we only take neutrino Yukawa couplings 
and the top Yukawa coupling into account. For the neutrino Yukawa coupling matrix 

3 Note that this statement is only true as long as there is no flavour symmetry dictating the form of 
the charged lepton Yukawa coupling Y e . 
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h v , this implies a system of nine differential equations, one for each element of /i„0 
As far as the input values are concerned, we choose them such that we start from a 
bimaximal mixing scheme [64] at the GUT scale. Although this mixing pattern is simply 
our assumption for the starting values at the input scale, it is nevertheless a well-motivated 
choice, as there are many flavour symmetries that tend to predict bimaximal mixing (see, 
e.g., Refs. |S51 E3 E3 EE1 E3 EZD1 EH)- At the GUT scale, this means that we choose our 
Yukawa coupling matrix to be given by 

/C put (M GUT ) = C^Mgdt) Vh-> ( 22 ) 

where we have taken the input values ^ iag (M G uT) = diag(0. 245, 0.42, 0.1), and £/bimax is 
the bimaximal mixing matrix, 



bimax 



tin 7J °N 

1 l l 

2 2 v/2 
V 2 ~2 72/ 



(23) 



Note that Eq. (1221) encodes the leptonic mixing correctly at the GUT scale, where we 
(arbitrarily) assume the charged lepton and heavy neutrino mass matrices to be diagonal. 
However, the running will in general introduce non-diagonal elements to all these matrices, 
which means that we have to make use of Eqs. ( 11 3 aft to 1 11 3d}) at lower scales. In particular, 
although it is perfectly justified to neglect the charged lepton Yukawa couplings in the 
FT at high scales, this is not valid anymore when we solve the equations for the effective 
operators. The reason for this is that there is a direct competition between A5 and Y e , cf. 
Eq. (I21bj) . Hence, in the effective theories, we also input the charged lepton masses in the 
RGE for k^ 22 \ Accordingly, we have taken into account the running of Y e in ET3, ET2, 
and ET1. However, as this running is not expected to be very strong, we have simply 
used the low-energy experimental values of the charged lepton masses as input at the high 
scale, and indeed their values turned out to hardly change in the evolution. 

As illustration for the running of the neutrino couplings, we have depicted the evolution 
of the diagonal entries of h u (not to be confused with /i^ iag ) in Fig. [5j We could also 
have presented a plot for the non-diagonal entries, but we prefer to focus on the actual 
masses and mixing angles, to be presented below, as these quantities are much easier to 
interpret physically. Note that, since some of the Yukawa couplings fail to be defined after 
certain heavy neutrinos have been integrated out (cf. Sec. EJ), their extrapolated evolution 



4 Note that we take the couplings to be real, for simplicity. In the general case, however, one would 
instead have nine complex equations, i.e., eighteen real ones, minus the number of phases that can be 
absorbed by redefinitions of the fermion fields. 
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is depicted by dashed lines in Fig. [5j As we can see from the figure, their running is 
relatively moderate. However, we have to take into account that all nine entries of h v will 
run, and that the structure of the light neutrino mass matrix, Eq. ( II Op . is also strongly 
influenced by the running of the heavy neutrino and scalar parameters. 

The RGEs for the scalar couplings, Eqs. ( I19a[) to fl!9ej) and (|20aj) to (120bl) . make use 
of those for gt, h u , and Y e . All these equations are interdependent, and therefore we have 
to solve the whole system at once. The initial values we choose are: 

Ai, 3l 4, 5 (M GUT ) = 10" 9 and A 2 (M GUT ) = 0.08, (24) 

where Mgut again denotes the GUT scale. These initial values were chosen such that 
three requirements are fulfilled: 

1. The correct relation between mi and the coupling Ai (obtained after EWSB) holds 
at the electroweak scale: v 2 = — 

2. Vacuum stability is ensured. 

3. The masses of the particles originating from the inert doublet remain consistent 
with collider constraints (i.e., m v ± < 200 GeV [63J). 

Although it would seem natural to choose the same input value for X 2 as for all the other 
couplings there arises a problem, as already mentioned in Sec. 14.41 Since A 2 decreases 
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Figure 6: The running of the scalar parameters. Note that m\ has to be negative, in 
order to guarantee a VEV for the SM-like Higgs, cf. Sec. 12.11 

with decreasing energy, we have to choose its input value large enough not to have A2 < 
at some point, i.e., not to break vacuum stability. This behaviour stems from the form 
of Eq. (119bj) . whose right-hand side has a strong tendency to be positive, as can be seen 
easily by realizing that the term 4\ 2 T U will be dominant as long as we have sizable Yukawa 
couplings h v which are nevertheless in the perturbative range. However, one might be able 
to improve on this point by studying a more general scenario. Still, the parameter choice 
in Eq. ( I24"j) is very interesting because, as expected, all couplings receive strong corrections 
with the exception of A5, which remains of O(l0~ 9 ). This behaviour, together with the 
(moderate) running of the mass parameters m 12 , is presented in Fig. El In addition to 
that, we have plotted the running of A5, drawn to a larger scale, in Fig. [7J Indeed, A5 
does run, but its running is weak due to the fact that A5 only receives multiplicative 
corrections, while the other couplings all have additive contributions. 

5.2 Masses and mixings 

Let us now turn to masses and mixing angles. First, we want discuss the evolution of the 
physical scalar masses, depicted in Fig. [HI This plot is obtained by simply inserting the 
runnings of all scalar parameters, cf. Fig. [61 into the explicit expressions for the masses, 
Eqs. Q). Due to the absence of a VEV for 77, the coupling A2 does not contribute to 
the scalar masses directly, which is why indeed the running of the inert scalar masses is 
mainly driven by the additive corrections to A3. The splitting between the electrically 
charged and the electrically neutral components of 77 arises due to the evolution of A4, 
while the splitting between the real and imaginary parts of 77 is only proportional to the 
small coupling A5, and hence it is not visible in the plot. Nevertheless, the smallness of 
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the mass different between 77^ and rfj is decisive for the correct neutrino masses, and a 
small splitting is enforced by the approximate lepton number conservation in the model. 

Note that we predict a relatively large value for the SM-like Higgs mass at low energies, 
rrih m 140 GeV. This is mainly due to the strong corrections to the parameter Ai, cf. 
right panel of Fig. [6j which arise from a complicated interplay of all other parameters. 
However, we are still very close to the range that is indicated by the newest LHC data 
from ATLAS [T2"| 175] and CMS [7U [73], so that we can hope to obtain similar (or even 
better) predictions by taking into account, e.g., some more Yukawa couplings. 

Next we turn to the heavy neutrino masses, i.e., the solutions of Eq. (]18]1 . The input 
values at the GUT were chosen such that we get consistent light neutrino parameters. 
The configuration we use is: 

(Mi, M 2 , M 3 ) = (7, 50, 100) TeV, (25) 

at the scale \x = Mgut- The corresponding evolution is depicted in the left panel of 
Fig. Again, these masses are not really defined in the model below the scales where 
the corresponding heavy neutrinos have been integrated out, and hence again the use 
of dashed lines in the plot. However, of course these particles still "exist", in the sense 
that they have visible contributions to low-energy physics, in particular the existence of 
non-zero light neutrino masses. It is only that these heavy neutrinos are kinematically 
decoupled. Note that the corrections to these masses actually seem quite moderate in 
the plot. One has to keep in mind, though, that the absolute corrections to the heavy 
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neutrino masses are indeed quite large (several hundreds of GeV), and it is just the high 
scale of their masses and the multiplicative nature of their corrections, cf. Eq. (fl8]) . which 
render their relative correction to be small. 

Finally, we want to discuss the light neutrino sector, which is probably the most 
important and interesting result in our simplified scenario. The evolution of the neutrino 
mass matrix can be determined by inserting the results of the relevant RGEs into Eq. ffTO]) . 
Hence, we actually do not choose directly any input values for the light neutrino masses 
at the GUT-scale, but we rather use the input values of the neutrino Yukawa and scalar 
couplings, as well as the ones of the heavy neutrino masses, which are all explicitly specified 
at the input scale. We can then determine the evolution of the neutrino mass matrix 
A4 u (fi), as a function of the energy scale /i. By calculating the eigenvalues of A4 U , we can 
determine the evolution of the light neutrino mass eigenvalues, depicted on the right panel 
of Fig. Note that we have chosen to work with a normal neutrino mass ordering, as 
explained above. We can see from the plots that the running of the light neutrino masses, 
and in particular the threshold effects, are relatively strong. This was to be expected, 
as there are many couplings involved in the physical neutrino mass matrix whose effects 
could potentially enhance each other. By taking a closer look at the plot, we can also 
recognize the seesaw structure: N 3 is the heavy neutrino with the largest mass, and it 
gets integrated out at the boundary between FT and ET3. Since, however, the inverse 
of the heavy neutrino matrix leads, by the seesaw structure, to the strongest suppression 
of the light neutrino masses, one can see the threshold effect for the lightest neutrino 
mass rrii at exactly that scale, while the other masses are hardly influenced. Similarly, 
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Figure 9: The running of the heavy and light neutrino masses. 



integrating out N2 results into a strong threshold effect for m2, and integrating out Ni 
strongly influences m^. 

Our goal is to reproduce the neutrino mass eigenvalues at low energies. However, as 
explained in the introduction, currently we still have no knowledge of the absolute neutrino 
mass scale, apart from upper bounds. Hence, we compare our low-energy predictions to 
the experimental values by "fixing" the light neutrino mass scale by the value obtained 
for mi at the Z-pole. We then take the values measured for the neutrino mass-square 
differences [TJ |2] to determine the experimentally preferred 3a ranges for m 2i3 , which are 
indicated by the black boxes in Fig. |9l right panel. Note that the experimental ranges 
are very narrow, and we only hit the corresponding range with rri2, while our prediction 
for ?ri3 is close to but still below the experimentally preferred range. On the other hand, 
as already explained, we are looking at a simplified scenario here, so we cannot expect to 
have a perfect match to the experimental values. Furthermore, we often measure the mass 
square differences at even lower energies (e.g. (9(MeV) for reactor neutrinos or 0(GeV) for 
atmospheric neutrinos [TB]), and our figure shows considerable running even below that 
scale. Apparently, it would be worth to perform a detailed investigation of this region in 
the Ma-model. 

Finally, the evolution of the three neutrino mixing angles is depicted in Fig. [TUJ It 
is this picture that illuminates the potential of our considerations most clearly: In the 
figure, the experimental 3a ranges for the mixing angles #12 (black) and #23 (red) [HE], 
as well as the newly measured 6> 13 (blue) [3] are indicated by the light colored boxes at 
the low-energy scale. This plot exhibits a remarkable behaviour of the Ma-model: We 
have started off with a bimaximal mixing pattern at the GUT scale, which could easily be 
motivated by a flavour symmetry!^] However, this mixing pattern is not compatible with 

Note that the simplest flavour symmetries have a certain preference for predicting maximal and zero 
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Figure 10: Evolution of the leptonic mixing angles. 



low-energy data. But what happens in our setting is that the evolution of the leptonic 
mass matrices is such that 6 ) 12 decreases and #i 3 increases with decreasing energy, while 6*23 
remains practically constant. Hence, we have a remarkably good match of our predictions 
with the low-energy experimental ranges. In particular, the corrections have a tendency 
to lead to a relatively large value for #13, which is in perfect agreement with the very 
recent measurements by Daya Bay and RENO. In that respect, our result is timely and 
suggests a more detailed study of the evolution of the parameters in the Ma-model, which 
we leave for future work. 

Note that, of course, our low-energy predictions strongly depend on the input values 
chosen at the GUT-scale. Had we chosen other parameters, we would have obtained 
different predictions. It had not been guaranteed that there is a certain tendency in the 
Ma-model in our setting not only to predict the correct mixing angles but also to be in 
agreement with the ranges or bounds for many other observables, but it actually turns 
out to be the case, at least in our example scenario. This is non-trivial, and it was not a 
priori clear that this would happen. 

Furthermore, there is an interesting test of the high-energy sector: One could apply 
different flavour symmetries to predict a set of mixing patterns at the GUT scale (or at 

mixing. 
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some other high scale) and evolve the neutrino observables down to low energies. Then, 
depending on the high-energy prediction, certain models would be in agreement with 
low-energy phenomenology, while others would violate certain conditions. In addition, in 
this setting, there are also alternative observables as, e.g., the inert scalar masses. This 
could even be extended to include, e.g., heavy neutrino or inert doublet Dark Matter or 
LFV processes. Hence, one could use the Ma-model as a tool to tell apart successful from 
failing symmetries: If a certain flavour symmetry predicts a high-energy pattern that does 
not yield low-energy agreement, we can discard it in the context of the Ma- model. This 
is particularly interesting since even relatively simple flavour symmetries, which predict 
a simple mixing pattern at high energies, could still be compatible with low-energy data. 
There is a rich variety of similarly interesting applications of our RGEs which could be 
studied in the future. 

6 Conclusions and outlook 

In this paper, we have investigated the 1-loop renormalization of Ma's scotogenic model, 
which is one of the easiest extensions of the Standard Model leading to a non-zero neutrino 
mass at 1-loop level. After having introduced the model itself, we have discussed how the 
different effective theories arise when integrating out the heavy neutrino mass eigenstates 
one by one. This results into three different effective theories at energies below the 
full theory. When calculating the renormalization, we need to match all operators that 
contribute to the neutrino mass at 1-loop level whenever a "new" effective theory arises. 
We have calculated the renormalization group equations for the full theory as well as for 
all three effective theories, which led to a very general expression for the light neutrino 
mass matrix in terms of various couplings present in the model. Furthermore, we have 
performed a numerical study of a simplified scenario, in order to give an illustration 
of the renormalization effects that arise in the model. As to be expected for a loop 
mass, the running and in particular the threshold effects are relatively strong: While the 
leptonic mixing angles exhibit considerable running at high energies, the running of the 
light neutrino mass eigenvalues originates mainly from threshold effects. We have found 
a particularly nice setting that suggests interesting connections to flavour symmetries: 
We start off with bimaximal mixing pattern at the grand unification scale, which evolves 
towards low energies into a pattern that matches very well the experimentally allowed 
ranges for the mixing angles. In addition, also the neutrino mass eigenvalues and the 
Higgs masses are at least close to the experimentally interesting regions. 

As we have seen, the renormalization of radiative neutrino mass models can involve a 
very rich and interesting phenomenology. While we have studied the case of the scotogenic 
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model, it could be similarly beneficial to compute renormalization group equations for 
other models with a loop-mass for neutrinos, such as the Zee-Babu model. In all these 
models, detailed phenomenological studies could be performed, in order to test the validity 
of the models thoroughly. It could be that the running is, e.g., too strong in certain cases, 
thereby tending to disagree with certain experimental bounds, or it could be that the Higgs 
sector in some of those models have severe problems with vacuum stability. In terms of 
neutrino physics, and in particular in the light of the recent measurements of $13, maybe 
one of the most interesting points to study is the combination of flavour symmetries with 
the known radiative neutrino mass models. We hope that many interesting such studies 
will appear in the future. 
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